Spectroscopic fingerprints of a surface Mott-Hubbard insulator: the case of SiC(OOOl) 
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We discuss the spectroscopic fingerprints that a surface Mott-Hubbard insulator should show 
at the intra-atomic level. The test case considered is that of the Si-terminated SiC(0001)\/3 x v3 
surface, which is known experimentally to be insulating. We argue that, due to the Mott-Hubbard 
' phenomenon, spin unpaired electrons in the Si adatom dangling bonds are expected to give rise to a 

Si-2p core level spectrum with a characteristic three-peaked structure, as seen experimentally. This 
structure results from the joint effect of intra-atomic exchange, spatial anisotropy, and spin-orbit 
coupling. Auger intensities are also discussed. 
Keywords: Semi-empirical models and model calculations; Photoemission; Auger electron spectroscopy; Magnetic 
phenomena; Surface electronic phenomena; Silicon carbide; Insulating surfaces. 
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The fractional adlayer structures on semiconductors show rich phase diagrams with potential instabilities of charge 
density wave (CDW) and spin density wave (SDW) type at low temperatures, as well as Mott insulating phases [Q. 
The (0001)\/3 x \/3 surface of hexagonal SiC, (as well as the closely analogous (lll)-\/3 x V3 of cubic SiC) expected 
from standard LDA calculations to be a 2D metal with a half-filled narrow band of surface states in the bulk energy 
gap - is instead experimentally proven to be an insulator with a rather large (2 eV) band gap Q,^]. It has 
been suggested that this system is a Mott insulator ||] due to the large value of ratio U/W, where U is the Coulomb 
interaction parameter, of the order of several eV, and W is the surface bandwidth, calculated to be about 0.35 eV 
gj. Very recent STM data further confirm this picture 
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A recent LSDA+U calculation pi] predicts an antifcrromagnetic insulating state with 3x3 magnetic periodicity and 
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a indirect gap of about 1.5eV, in reasonable agreement with the experimental value J^g. An exchange coupling J 
between neighboring Si adatoms of about 30K is also calculated. The LSDA+U magnetic state, with spin-collinearity 
enforced in the calculation, turns out to have a uniform magnetization m z — 1/3: one spin-down and two spin- 
up adatom surface bands are completely filled in the 3x3 Brillouin zone. Crudely speaking, this corresponds to 
two adatom dangling-bond orbitals having spin-up electrons and one having spin-down electron in each triangular 
__ ■ plaquette. When spin non-collinearity is allowed one finds, e.g., from Hartree-Fock Hubbard model calculations, ||] 
\& \ that the state with spins lying on a plane - in a 120° three-sublattice Neel structure, typical of the triangular lattice 
Heisenberg antiferromagnet - has slightly better energy than the collinear one, but the resulting energy gain is so 
small, about 0.5 meV/adatom, that unrealistically precise ab-initio calculations would be necessary to decide which 
of the two is the actual ground state. Moreover, the bands for the two solutions are rather similar, apart from some 
| extra splittings introduced by the collinear magnetic solution 

In reality, both the LSDA+U and the Hartree-Fock Hubbard calculations, being mean-field approaches are at best 
only cartoons for the actual Mott-Hubbard state, which is insulating not because it is magnetic, but rather due to 
strong correlations inhibiting electron hopping. Once this is appreciated, one realizes that - even in the absence of 
any true magnetic long-range order - finite magnetic moments on the Si adatoms should be experimentally detectable 
up to high temperatures, so long as there is an insulating gap. The low-temperature low-energy collective physics 
of these unpaired adatom spins is instead governed by a much smaller energy scale, the exchange coupling J. Is it 
O ■ conceivable that a genuine low-temperature spin long-range order (LRO) would be detectable in this quasi-2D system? 
It would seem unlikely, since: 1) The estimates for J already suggest that LRO, if present, must be there only at very 
• rH . low temperature. Moreover, thermal fluctuations might destroy LRO altogether, if the continuous spin symmetry is 
not broken by anisotropy effects due to spin-orbit coupling; 2) It is at present unclear which spin ordering will be 
stabilized by the spin-orbit coupling; 3) If spins lie on a plane, with zero net magnetization, they are presumably 
hard to detect with LRO sensitive techniques. Therefore, without ruling out totally the possibility of low-temperature 
LRO, we believe that the SiC(0001)\/3 x -\/3 surface will most likely be in an overall paramagnetic Mott insulating 
state, in spite of the existence of an on-site moment, at least down to liquid nitrogen temperature. 

It is the aim of the present paper to discuss in some detail what kind of experimental signatures one should look 
for in a surface Mott-Hubbard insulator, apart from a gap in the photoemission spectrum. As mentioned previously, 
methods which require magnetic LRO are likely not viable. The situation could be much more promising at the 
intra-atomic level. We will argue below that the very fact that the system is in a Mott-Hubbard insulating state, 
independently of the magnetic LRO issue, brings about unpaired spins in the adatom dangling bonds, which in turn 
determine exchange splittings in the adatom core levels which should be large and detectable. These splittings may 
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in fact have been already observed for SiC(OOOl) in the surface component of the Si-2p core level photoemission 
spectrum, although attributed so far to some unspecified charge inequivalence of two types of Si adatoms. 

To make things more concrete, let us assume that, due to the Mott-Hubbard phenomenon, the 3p z -like dangling 
bond orbital of the Si adatom in the V3 x \fi SiC(OOOl) is singly occupied by an unpaired electron, say of spin up. 
In reality, the Wannier function for the surface state is not a pure 3p z , since it extends itself down into the first SiC 
bilayer the weight a of the Si 3p z orbital in the Wannier function is about 0.5, the remaining lobe having a 

negligeable overlap with the Si-adatom core The 3s and 3p XiV orbitals of the Si-adatom are involved in the 

bonding with first-layer C atoms, and contain spin-paired electrons. Therefore, a Si-adatom has two distinct sources 
of anisotropy: i) the first comes about because the occupation a of the 3p z orbital is not equal to the occupation, call 
it f3, of the 3p x or 3p y orbitals; ii) the second, and most important, source of anisotropy is due to the fact that the 
electron in the 3p z dangling bond is spin unpaired, as opposed to the standard singlet bonding configuration in 3p x . y . 
These anisotropics break the core level symmetry between 2p z and 2p xy by an amount which we will argue to be not 
negligeable compared to spin-orbit. We will now estimate the expected 2p core level splittings, with and without a 
spin-orbit coupling, and show that while spatial anisotropy of type i) by itself would result in a very small effect, the 
spin induced anisotropy ii) has a large influence on the spectrum. 

The calculation we perform is very simple, and is based on the Si 2p — 3p Slater exchange and Coulomb integrals, 
obtained from pure atomic Hartree-Fock calculations. Assume the Si adatom valence electrons to be distributed with 
weights a in 3p z j, and (3/2 in each 3p x>y | and J.. Let us consider first the situation in absence of spin-orbit coupling, 
Aso = 0. The Hartree-Fock energy of the 2p z j states can be written as 

(3 

(■2p z \ = a [ K 2p z ,3p z - J2p z ,3p z ] + 2— ^i^p^Sps - J2p z ,3pJ + ■ ■ ■ , (1) 

where the dots stand for one-body terms and two-body contributions due to other spin-paired occupied or- 
bitals which will exactly cancel out in all the energy differences which we will calculate. Here K 2pz3p ^ x - ) = 
(2p z ,2,p z(x) \e 2 I\yi - r 2 \\2p z ,3p z[x) ), and J 2pz ^ Pz{x) = (2p z ,3p z{x) \e 2 /\y x - v 2 \\3p z(x) ,2p z ) , are the relevant 2p - 3p 
Coulomb and exchange Slater integrals. The factor 2 in front of the (3/2 term is due to equal contributions coming 
from 3p y and 3p x . The factor 2 in front of K 2pz ^ P:c keeps into account the direct contribution due to both ] and [ 
3Px,y electrons. (The corresponding exchange term has no factor 2 because only 3p x . y j electron contribute to the 
exchange gain of the 2p z ] level.) A similar calculation for £2 Px j = tip v x gives 

(3 (3 

£2p x T = a[^2p x ,3p 2 - ^2p x ,3pJ + -^[2-?Gp x ,3p x - ^2p x ,3pJ + ~Z [2-ftT 2 p x ,3p H - <^2p x ,3p B ] + ■ ■ ■ , (2) 

which can be reexpressed, due to rotational symmetry, in terms of the previously introduced K 2pz}3pz ( x } and J 2pz ^, Pz ( x )- 
We now estimate the Slater integrals assuming pure atomic hydrogen-like 2p — 3p orbitals (the effect of intra-atomic 
hybridization in a full atomic HF calculation would change a bit the numbers: these effects will be estimated later on). 
Numerically, the two different Coulomb integrals appearing in Eqs. |l|-|| are large and comparable (K 2pz ^ Pz ss 30.5 
eV, K 2pz ^ P:L ~ 28.2 eV), whereas the 2p z — 3p z exchange integral is smaller (J 2pz ^ Pz ~ 4.3 eV), but one order of 
magnitude larger than the 2p z — 3p x exchange integral (J 2pzt 3 P:L ~ 0.5 eV). The large difference between J 2pz ^ Pz and 
J 2 p Zi 3 Px immediately implies that the exchange splitting between 2p z f and 2p z J. is rather large 

£2p 2 j. — e 2 Pi T ~ o:J 2pz ^ Pz » 2eV 
whereas the correponding exchange splitting between 2p x j y \ and 2p x j y \ is one order of magnitude smaller: 

£ 2p x/y l - e 2 Px/M T ~ otJ 2pz , Zpx « 0.2eV 
We can also calculate, using Eqs. |l[|^, the splitting A = £2p x/M T — e 2p^T between the 2p z f and the 2p x / y f levels, 

A = a[(J 2pzi3pz - J 2 p i: 3pJ - {K 2pz ^ Pz - ^2p,,3p x )] + ^[2(^2p^3p, - ^2p,,3p x ) ~ {J2p z ,3p z ~ ^2p z: 3p x )] ■ (3) 

Notice that the structure of the a contribution in Eq. [| differs from that of the (3 term because the 3p z is assumed 
to be spin unpaired (anisotropy ii), above): If we were to assume that the occupancy a was due a/2 3p z | and a/2 
3p z |, the a term would be given by {a/2)[{J 2pz ^ Pz - J 2pz ,3 P!C ) - 2(K 2pz ^ Pz - K 2pz ^ P:c )] and would exactly cancel 
off the f3 term when a = (3 (bulk case). Plugging in the numbers for the Coulomb and exchange integrals we get 
[2(K 2pzi3pz - K 2pz ^ Px ) - (J 2pzi 3 Pz - ^2p,,3pJ] ~ [2(2.3eV) - 3.8eV] = O.GeV, and [(J 2 p z ,3 P , - J2 Pz ,3 P J ~ (K 2pz ,3 Pz - 
#2p„3pJ] « [3.8eV - 2.3eV] = 1.5eV, which provide the desidered estimate for A » a[1.5eV] + (P/2)[0.6eV], i.e., 
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A w leV for a w 0.5 and /3 w 1. The effect of intra-atomic orbital hybridization which modifies the atomic Hartree- 
Fock orbitals with respect to pure hydrogen-like 2p — 3p atomic orbitals, is estimated, from a full Hartree-Fock Si + 
atomic calculation, to further reduce A and the 2p z / x exchange splittings by a factor of about 2, so that we arrive at 
the crude estimate A w 0.5eV, e 2pz i — ^2p,\ ~ leV - , and e 2pm i — ^2p x f ~ Q.leV. The overall picture for the Si-adatom 
2p spectrum in absence of spin-orbit, but with all the previously described anisotropy effects included, is summarized 
in Fig. 0(b). When we include the spin-orbit interaction AsoL • S, the final core hole levels are obtained as in Fig. 
0. For Aso = 0.4 eV (the value extracted from the bulk 2p J=1 / 2 — 2p J=3 / 2 experimental core level splitting) we 
find a roughly three-peaked structure, the broad central peak four times as strong as each side peak. This offers an 
alternative explanation of the experimental lineshape pi (see Fig. ||) in terms of a single exchange- split multiplet, 
rather than two chemically inequivalent sites Si and S 2 p0[ whose existence is otherwise not supported. Incidentally, a 
similar three-peaked structure for the Si-2p core level spectrum has been experimentally observed for a sub-monolayer 
Si segregated at the surface of iron |IlJ . 

Additional direct experimental evidence for the Mott-Hubbard state could be obtained by a careful study of adatom 
(3p 2 , 3p z ) Auger spectral intensities, which should be easily singled out owing to the large gaps. The probability of 
double occupancy of the adatom dangling bond orbital should be almost completely suppressed, dropping from the 
band value of 1/4, to a value of order t/U - where t is the hopping matrix element and U is the effective on-site 
Coulomb repulsion - which is one order of magnitude smaller. Even considering that only half of the orbital is adatom 
3p z , this surface should show a {2>p z , 3pz) Auger intensity which, by comparison with the remaining (3p, 3p) values, 
is anomalously small, as a proof of its Mott-Hubbard state. 

Summarizing, we have argued that the Mott-Hubbard insulating surface state in SiC(OOOl) should leave clear 
and detectable intra-atomic fingerprints. The Si 2p core hole will exhibit a 6-fold multiplet resulting from the joint 
effect of intra-atomic exchange, asymmetry, and spin-orbit, with a characteristic three-peaked structure. Auger 
intensities should also be affected. More experimental and theoretical effort is clearly called for to check these strong 
correlations and related magnetic effects, possibly the first of this magnitude to be suggested for an sp-bonded, valence 
semiconductor surface. 
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FIG. 1. Schematic 2p core level energy scheme for a bulk Si atom (a), 
Estimates for the splittings in (b), discussed in the text, are A ~ 0.5eV, J zz : 



and for a spin- unpaired surface Si adatom (b). 
s leV, and J zx w O.leV. 




FIG. 2. Si 2p core levels calculated as a function of spin-orbit coupling, in comparison with photoemission data and original 
fitting in terms of bulk (Bl) and two surface sites SI and S2 by Johansson et al. (lo| . Upon inclusion of intra-atomic exchange 
splitting the whole surface contribution S1+S2 can be explained as due to a single site. This figure is taken from Ref. M. 
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